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Abstract. The Unified Radio and Plasma Wave (URAP) experiment on Ulysses has 
provided unique high latitude measurements of Jovian hectometric radiation (HOM) 
during its encounter with Jupiter in February 1992. URAP was the first radio instrument 
in the Jovian environment with radio direction-finding capability, which was previously 
used to determine the HOM source locations in the Jovian magnetosphere. These initial 
source location determinations were based on several assumptions, including the neglect 
of refractive effects, which may be tested. We have, for the first time, combined the 
measured incident ray-direction at the spacecraft with a model magnetosphere to directly 
trace the rays back to the HOM source. We concentrate on the observations of HOM 
from high northern latitudes when Ulysses was at distances <15 R r The three- 
dimensional ray-tracing calculations presented here indicate that the HOM sources 
probably lie on L shells in the range 3sL<7 (tilted dipole magnetic field model) 
consistent with previous determinations that ignored the effects of refraction. The ray- 
tracing results, however, indicate that wave refraction due to the Io torus and the 
magnetic field can significantly influence the precise source location. We show that 
constraints on the locations imposed by the gyroemission mechanism suggest that the Io 
torus density may have experienced temporal and/or spatial fluctuations during the Ulysses 
observations of HOM. Finally, in the cold plasma approximation we demonstrate that 
even if the emission were nearly linearly polarized near the source region, almost circular 
polarization will be observed at Ulysses, in agreement with observations. 


1. Introduction 

Intense hectometric radio emission (HOM) is observed at 
Jupiter and lies in the frequency range from perhaps 100 kHz 
to over 2 MHz [cf. Garret al. , 1983]. At the higher frequencies 
it is difficult to distinguish from decametric (DAM) radio emis- 
sion, which has been well studied. In fact, recent observations 
of Jovian emissions made by the Wind satellite suggest that 
perhaps much of HOM is a low-frequency extension of DAM 
(M. L. Kaiser, personal communication, 1996). HOM is most 
probably RX mode emission generated by the cyclotron maser 
instability. 

The source location of HOM has been controversial for 
many years. Calvert [1983] analyzed the Faraday rotation of the 
HOM radio emission observed by Voyager 1 near closest ap- 
proach and concluded that this radio source was located near 
the gyrosurface and along L = 6 (the Io flux tube). Ladreiter 
and Leblanc [1989, 1990a, b| have reported ray-tracing studies 
of HOM observed by the Voyager flybys of Jupiter. Ladreiter 
and Leblanc [1990a] argued that ray-tracing calculations, per- 
formed with large wave normal angles (70°— 90°) as suggested 
by the cyclotron maser instability mechanism, indicate emis- 
sion sources located along field lines in the range of L shells 
15 < I. < 20 (Jovian auroral zone). From these source 
locations, the ray paths to the Voyager spacecraft (located near 
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the Jovian equatorial plane) will be significantly refracted by 
the Io plasma torus, as a function of frequency, resulting in an 
emission shadow zone that may extend radially 20 R, or more. 
Since Ladreiter and Leblanc had no a priori knowledge of the 
direction of propagation of the incident radiation, they consid- 
ered only rays in the meridian plane propagating parallel to the 
magnetic equator. For ray propagation directed out of the 
meridian plane, the source locations may be significantly altered. 

Reiner et al. [1993a, b] have reported the direct determina- 
tion of HOM source locations based on direction of arrival 
measurements made by Ulysses during the Jupiter flyby. In 
contrast to Ladreiter and Leblanc [1989, 1990a, b], knowledge 
of the incident direction of the radiation provides a direct, 
clearly defined constraint on the possible locations of the 
HOM sources in the Jovian magnetosphere, without the ne- 
cessity of introducing constraints from the emission mechanism. 

The Ulysses HOM direction-finding results were used to 
determine HOM source locations for two different observa- 
tional geometries. Reiner et al. [1993b] reported the determi- 
nation of HOM sources detected when Ulysses was near the 
equatorial plane at distances r > 40 R r These results indi- 
cated both northern and southern hemisphere HOM sources 
all located on low L shells, 4 < L < 6. Reiner et al., [1993a] 
reported observations of HOM made while Ulysses was lo- 
cated at r < 14 Rj and at Jovigraphic latitudes >30°. in this 
work the HOM sources were again found to be located at low 
L shells along field lines that could pass through the Io torus. 
In subsequent studies, Ladreiter et al. [1994], using a slightly 
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Figure 1. A projection of the Ulysses spacecraft trajectory 
superimposed on contours of the lo plasma torus density as 
determined from the model by Divine and Garrett [1983]. The 
locations of events 1 and 2 are indicated. 


different method of analysis, found auroral region sources for 
HOM in the range 7 < L < 11. For ail of these analyses, 
straight-line propagation of the radiation was assumed, and 
possible refractive effects were ignored. 

In this paper we have combined the measured incoming 
wave direction with ray tracing in a model Jovian magneto- 
sphere. The model magnetosphere combines the densities [Di- 
vine and Garrett , 1983] with the 06 magnetic field model [Con- 
nerney , 1993]. The location of the radio source is determined by 
tracing the rays from the spacecraft back along the direction of 
the incident radiation to the RX cutoff'. The results of this 
paper are consistent w'ith HOM sources located at low L shells 
but also indicate that the Divine-Garrett time-independent 
model may not be completely adequate to account for the 
HOM source directions observed during the February 1992 
Ulysses encounter. 

2. Model 

The ray-tracing code is based on the Haselgrove [1955] equa- 
tions and has been discussed in the past [cf. Menietti et al ., 
1987]. For the Jovian magnetosphere we have incorporated the 
analytical plasma density model of Divine and Garrett [1983], 
and the 06 magnetic field model of Connerney [ 1993]. The 06 
magnetic field model is a fit to the Pioneer 1 1 and Voyager 
observations using a spherical harmonic expansion and is also 
a static model. This model is an improvement over the 04 model 
[Acuna and Ness , 1976] as discussed by Connerney [1993]. 

The Divine-Garrett plasma torus model is a three- 
dimensional (3-D), static, analytical, and empirical model with 
azimuthal symmetry and is based on the Pioneer and Voyager 
flyby data. It incorporates the Io plasma torus densities first 
observed in detail by the Voyager spacecraft [Bagenal et aL , 
1985; Bagenal and Sullivan , 1981]. It should be noted that no 
direct measurements of plasma have been made near the 
planet (r < 4 R f ) , and the Divine-Garrett model in this region 
is an extrapolation. The torus is believed to consist of separate 
cold and warm regions, and the boundary between these re- 
gions extends to perhaps 3 R f along magnetic field lines. This 


empirical model was used because it lends itself well to ray- 
tracing calculations and agrees reasonably well (within a factor 
of ~2) to the most recently published tabular model [ Bagenal , 
1994]. In Figure 1 wc display contours of the plasma density of 
the Io torus determined from the Divine and Garrett [1983] 
model (modified version of their Figure 1 1). 

It is becoming increasingly clear, however, that torus models 
based on the Voyager flyby observations may not adequately 
account for spatial and temporal variations of plasma density. 
Desch et al. [1994], for instance, conclude that significant azi- 
muthal asymmetries (namely, 30-50% density depletions) are 
necessary to explain Ulysses radio observations in the high 
latitude Io plasma torus. Schneider and Spencer [1991] have 
obtained ground-based observations indicating strong longitu- 
dinal and temporal dependences associated with the Io torus, 
which no doubt influence radio emissions that arc refracted by 
the high density regions. 

The ray-tracing code requires as input the initial wavenor- 
mal angle (i j/ () ) and the azimuthal angle, j3 0 , of the k vector 
about the magnetic field direction. These angles were derived 
directly from the Ulysses observations at the spacecraft as 
opposed to the HOM source point. The measured direction 
angles at Ulysses, which define the direction of the ray wave 
vector k at the spacecraft, were first transformed from the 
spacecraft frame lo Jovian System III coordinates. To deter- 
mine the wave normal angle of the ray from these direction 
angles, a transformation was made to a coordinate system 
directed along the B field vector direction at Ulysses. The wave 
direction angles at Ulysses are defined in this system where ^ 
is the angle between B and k (the wave normal angle) and is 
the azimuthal angle of the projection of k on the plane per- 
pendicular to B. 

The direction-finding capability makes use of a synthesis of 
signals from the spin plane (72.5 m tip-to-tip dipole) and spin 
axis (7.5 m monopole) antennas. These data provide the arrival 
direction and complete polarization state (four Stokes param- 
eters) of the radiation [Manning and Fainherg , 1980]. The two 
angles that identify the direction of the incoming HOM radi- 
ation were obtained from the analysis of the resulting spin 
modulated signals from the radio receivers observed for each 
12-s spacecraft spin period. The sampling rate is frequency 
dependent. Variations of the source intensity over the 12-s 
measurement cycle, as well as errors due to receiver noise and 
quantization errors in the telemetry, lead to a significant scat- 
ter in the individual angle determinations [e.g., Reiner et «/., 
1993a, Figure 2]. To reduce these errors, the source directions 
used in the ray tracing were obtained by averaging the individ- 
ual direction measurements over a period of about 30 min. 
This averaging introduces errors due to the motions of the 
spacecraft and the source in the time over which measure- 
ments were made, but these latter errors are estimated to be 
significantly smaller than those due to the above effects. Fi- 
nally, there may be systematic errors in these measured angles 
due to the calibration of the receiver/antenna system. To cal- 
ibrate, we used the observations of HOM in the same fre- 
quency range measured when Ulysses was far from Jupiter 
(>50 Rj) to determine the receiver gain and electrical tilt 
angles of the monopole spin axis antenna. We estimate that 
this should allow us to determine the source angles to ~ 0.5° in 
the spacecraft system. 

The rays al each frequency were launched from the location 
of Ulysses at the measured direction angles back toward Jupi- 
ter and followed at least until the ray coordinates reached a 
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Table I. Observed HOM Emission Angles and Deduced 
Source Locations 



Direction Angles 

Source Location 

Frequency 

kHz 

deg 

ft„ L Shell 

deg (Dipole) 

r(Rj) 

A, 

deg 

Am, 

deg 

540 

153.8 ± 2 

Event / 

262.0 ± 4 5.3 ± 0.2 

3.20 

39.0 

49 

940 

151.6 ± 2 

270.0 ±4 3.9 ±0.1 

2.62 

35.0 

34 

540 

153.4 ± 2 

Event 2 

217.1 ±4 7.2 ±0.2 

3.30 

47.4 

105 

740 

151.4 ± 2 

223.8 ±4 6.4 ± 0.3 

2.97 

47.1 

116 


No refraction [cf. Reiner et a /., 1993a]. 


point (“source point”) where the wave frequency was n ear the 
local RX cutoff frequency, / RX = fj 2 + V (/ p ) 2 + (fJ2) 2 , 
where f s is the gyrofrequency and / p is the plasma frequency. 
As a practical criterion, we considered the source position to 
be any point within 5% of the local RX cutoff frequency. 

3. Results 

We have concentrated on events 1 and 2 as presented by 
Reiner et al. [1993a]. For these events, Ulysses was located at a 
distance of about 12 Rj and 6 Rj from Jupiter, respectively, 
and at latitudes >30°. Event 1 occurred on February 8, 1992, 
from about 0100 to 0300 spacecraft event time (SCET) when 
Ulysses was at system III longitude, A m ~ 31°, and latitude, 
A m 31°, while event 2 occurred from about 1000 to 1130 
SCET while Ulysses was at A n[ — 41° and A m ~ 37°. Specif- 
ically we have determined the wave direction angles listed in 
Table 1 for the SCET times 0215 (event l)and 10:30 (Event 2). 
As indicated, these angles have an associated error that results 
primarily from uncertainties in fitting the data to the antenna 
equations as briefly discussed above and elsewhere [cf. Reiner 
et al. , 1993a]. Included in Table 1 are the source positions 
obtained assuming no refraction of the ray [cf. Reiner et al ., 
1993a] for comparison with the determinations in this paper. 
The positions of the Ulysses spacecraft at the times of events 1 
and 2 are superimposed on Figure 1. 

3.1. Event 1 

We launched a ray at the angles \p 0 and ft, given in Table 1 
for / = 540 kHz, and traced this ray through the model Jovian 
magnetosphere toward Jupiter. Assuming that HOM emission 
is generated in the RX mode, this ray must intercept / RX at 
some point along the ray path, if we have chosen the initial 
wave direction correctly and if the magnetospheric model is 
valid. In Figure 2 we display the ray path in three dimensions 
(Figure 2a) and in a two-dimensional projection into the x-y 
plane (Figure 2b). Included on the plot are a portion of the 
spacecraft trajectory, the magnetic field line loop (06 model) 
that passes through the spacecraft, and three line profiles rep- 
resenting planar cuts through the surface for / RX = 540 kHz. 
The star on the plot indicates the projection of the 540-kHz 
RX cutoff surface at a position along the z axis that corre- 
sponds to the z coordinate of the ray path at the point of 
closest approach to Jupiter. We see that as the ray propagates 
toward Jupiter, it is refracted out of the meridian plane and 
does not reach the RX cutoff surface for / = 540 kHz but is 
in fact about 0.25 Rj above this surface at the closest point 



Figure 2a. Three-dimensional view, using system III coor- 
dinates, of the ray path, spacecraft trajectory, 06 field line 
through the spacecraft location, and planar cuts through the 
540-kHz RX cutoff surface, for the observed initial wave di- 
rection of event 1 . 

(/rx = 412.5 kHz). Since this particular ray refracted away 
from the planet just above the RX cutoff surface, / RX = 540 
kHz, either the initial wave direction was in error or the mag- 
netospheric models for plasma density and magnetic field did 
not accurately represent the Jovian environment at the time of 
the Ulysses flyby. 

By slightly adjusting the launch angles at the spacecraft, 
however, we were able to produce a ray path that did intercept 
/rx (within 5%). This was accomplished in a systematic fashion 
by alternately incrementing ip then ft If ip was increased and 



the 06 magnetic field loop that intercepts the spacecraft, and 
three line profiles of points on the 540-kHz RX cutoff surface. 
The star locates the projection of the 540-kHz RX cutoff sur- 
face at a position along the z axis where the wave ray path is 
closest to Jupiter. 
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Table 2. Event 1 HOM Launch Angles and Calculated 
Source Locations for 06 Magnetic Field Model 


Direction 

Angles 

L 

Torus 

Source Location 

Ak. 

deg 

deg 

P* 

deg 

r, 

Kj 

A, 

deg 

An,. 

deg 




/ = 540 kHz 




152.8 

262.0 

4.0* 

yes 

3.12 

34.6 

72.0 

1.68 

153.8 

269.0 

5.5f 

yes 

3.22 

48.8 

48.3 

3.09 




f = 940 kHz 




149.9 

269.8 

2.3* 

yes 

2.31 

16.1 

38.8 

1.68 

153.9 

270.8 

6.5t 

yes 

2.78 

58.5 

35.0 

2.35 

*/ 

imir 


the ray was found to move away from / Rx , then ijj would be 
decreased. Then p would be incremented in the same fashion. 
The ray, with initial wave direction angles at Ulysses of <// (l = 
152.8°, and p () = 265.0° intercepted the “source point” where 
the wave frequency was approximately equal to/ RX atL = 4.0 
(tilted dipole model). The system III latitude and longitude of 
the source point is tabulated in row 1 of Table 2. Also included 
in Table 2 is the angular difference, Ak, between the wave 
direction measured at Ulysses (refer to Table 1) and that used 
in the ray tracing. These values are generally smaller than the 
maximum estimated RMS error (—4°) in measuring k (which 
depends on i// () and j3 0 listed in Table 1). The results in row 1 
of Table 2 are for a ray as close as possible to the measured 
direction (Table 1) that reached / RX . It was also the ray that 
produced the smallest value of L (L min ) at the source location. 

To gain some idea of the range of possible source locations, 
we next varied values of t//,, and /3 0 until w’e obtained a ray that 
intercepted the f RX surface at a point with the largest possible 
value of L. The wave direction angles at Ulysses were i// (l = 
153.8°, j3 (l = 269.0 and the resulting source locations are tab- 
ulated in row 2 of Table 2. In Figure 3 we show plots of the two 
rays identified in rows 1 and 2 of Table 2. The ray paths are 
shown in a three-dimensional projection (Figure 3a), and in 
the x-y Jovigraphic plane (Figure 3b) as a “top” view. The two 
rays (one for L min and one for L max ) propagate out of the 
meridian plane but intercept the RX cutoff surface relatively 
near the spacecraft longitude. 

Table 2 also gives the Jovian system III coordinates of the 
“source point” of the emission or position along the ray path 
where the wave frequency was within 5% of the local value of 
/rx- By comparing these with the values in Table 1, we see that 
the source position determined from the ray tracing can differ 
considerably from that determined assuming no refraction and 
a dipole field. The range of source latitudes on the / RX surface 
deduced from the ray tracing is from 34.6° to 48.8°. The source 
longitude ranged from 48.3° to 72.0°. These results indicate 
considerable uncertainty in the HOM source locations. The 
values of 39.0° and 49.0° from Table 1 lie within this range. 

Also contained in Table 2 are the results for the ray tracing 
at a frequency of 940 kHz. By comparing with Table 1 we see 
that the values of the source position for the minimum and 
maximum L value bracket the values of Table 1 for no refrac- 
tion; and the longitudes obtained by ray tracing arc also quite 
close to those shown in Table 1. This result is not surprising 
considering the smaller refractive effects expected at the higher 
frequency. 


3.2. Event 2 

Event 2 occurred near closest approach of Ulysses at a radial 
distance of ~6 R f . As seen in Figure 1, this event occurred 
when Ulysses was near the upper edge of the torus. For event 
2, the ray corresponding to the measured w'ave direction angles 
again failed to propagate to the RX cutoff surface. To intercept 
the f RX surface, it was necessary to considerably modify the 
values of i// 0 and p {) from their measured values. The results of 
the analysis for event 2 at / = 540 kHz and/ = 740 kHz are 
presented in Table 3. The source locations for event 2 are 
similar to those of event 1 but generally occur at smaller lon- 
gitudes than those determined assuming no-refractive effects 
(Table 1). We have included in row I of Table 3, as in Table 2, 
the results of the ray launched as close as possible to the 
measured direction that intercepted the RX cutoff surface. The 
source location for this ray (the point where f ~ / RX ) was 
located at the smallest value of L for all rays examined. In row 
2 we include the results for the ray launched in a direction that 




Figure 3. The format is the same as Figure 2, but now the 
rays launched in directions that yielded source positions at the 
smallest and largest possible L values are shown. The ray 
labeled L min is also the ray launched as near as possible to the 
observed ray direction. 
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yielded a source location at the largest possible value of L 
(/-max)- Both of these rays (rows I and 2 of Table 2) are plotted 
in Figure 4a (top view) and Figure 4b (3-D view). 

As seen in Table 3, the values of Ak are significantly larger 
than the estimated measurement uncertainty in the arrival 
directions of the radiation. To understand the origin of this 
result we analyzed the sensitivity of ray paths to modifications 
in the plasma density model. We have both increased and 
decreased the plasma torus densities by multiplying the values 
produced by the Divine and Garrett model by a factor ranging 
from 0 to 5. We found that we could obtain a smaller value of 
Ak by decreasing the torus density. In Table 3 we also present 
the ray-tracing results obtained in the absence of a plasma 
torus (including only a plasmaspheric density model). In row 3 
of Table 3 we indicate the results of the ray that produced not 
only the largest value of L but also the smallest value of Ak. 
The latitude of the source location of this ray lies intermediate 
between that of rows 1 and 2 in Table 3 (which were obtained 
including the Io torus), and the longitude is only ~ 4° larger. 
Even though the value of Ak is smaller than the results that 
include a torus ( rows 1 and 2 of Table 3), it is still considerably 
larger than the expected maximum observational error of 
about 4°. 

Another parameter in the problem is the Jovian magnetic 
field. In order to investigate the influence of modifications of 
the field model, we repeated many of the calculations at a wave 
frequency of / = 540 kHz using the 04 magnetic field model 
[Acuna and Ness, 1976]. For a detailed comparison of both 
models, see Conmrmy [1993]. In Table 4 we present the cal- 
culations of ray tracing for both events 1 and 2 in the same 
format as Table 3. The source locations deduced for the 04 
magnetic field model are qualitatively similar to those for the 
06 model. As in Tables 2 and 3, in row 1 of Table 4 we include 
the source location for the ray launched as close as possible to 
the measured direction, while row 2 contains the source loca- 
tion at the largest possible L value. For the 04 model the 
range of longitude values is somewhat smaller compared to the 
06 model, while the latitude range is larger extending to 57.5° 
at L m;ix . For event I the values of Ak were less than the 
expected measurement error. 

Table 4 indicates that while the 04 model produced a 
smaller range of longitudes for the source region than did the 
06 model, the values of Ak were still larger than the expected 


Table 3. Event 2 HOM Launch Angles and Calculated 
Source Locations for 06 Magnetic Field Model 


Direction 

Angles 


Source Location 


deg 

ft, 

deg 

L 

r. 

Torus Rj 

A, 

deg 

A,„. 

deg 

Ak, 

deg 

144.7 

249,5 

3.8* 

/ = 540 kHz 
ves 3.10 

31.3 

78.3 

18.5 

149.7 

271.5 

4.7i 

ves 3.13 

44.4 

42.9 

25.4 

146.7 

240.5 

4.9t 

no 3.25 

40.7 

r i 
x 

13.4 

143.8 

255.8 

3.2 

f - 740 kHz 
yes 2.73 

28.7 

72.4 

18.5 

150.8 

260.8 

5. Of 

yes 2.95 

46.7 

70.9 

17.6 

140.8 

248.8 

3.3 

no 2.76 

30.0 

77.3 

17.3 

144.8 

249.8 

3.94 

no 2.85 

36.9 

77.4 

15.1 


*L mil 




Figure 4. (a) Top and (b) three-dimensional view of the 

“L min " and “L Iliax ” rays for / = 540 kHz w hen the spacecraft 
was at the position of event 2. For this event there was a bigger 
range of possible source positions, and Ak was larger than the 
expected error. 

measurement error of —4°. Also contained in Table 4 are the 
results for the case of no torus density. The range of allowable 
ray paths have initial wave directions that lie closer to the 
measured values, but Ak is still larger than 4°. For comparison, 
we show in Figure 5 sample ray paths for Event 2 using both 
the 06 and the 04 magnetic field model. The rays were each 
launched in the same direction relative to the field line. The 
small displacement of the 04 model relative to the 06 no 
doubt accounts for the larger longitude values obtained in the 
former ray-tracing results. 

4. Polarization 

Recent observations of DAM [Dulk el al . , 1994, and refer- 
ences therein] indicate that this emission is often elliptically 
polarized and that the axial ratio may depend on the viewing 
angle [Melrose and Dulk , 1993]. Yet Reiner et al. [1993a] re- 
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Table 4. Modeled HOM Launch Angles and Source 
Locations for 04 Magnetic Field Model 


Direction 

Angles 

L 

Torus 

Source Location 

Ak, 

deg 

‘/'th 

deg 

Ah 

deg 

r, 

R, 

A, 

deg 

Am. 

deg 




Event 1: f 

= 540 kHz 



151.6 

259. 1 

3.2 

yes 

2.89 

24.9 

64.7 

3.72 

155.1 

263.5 

7.5 : 

* yes 

3.28 

57.5 

43.5 

4.01 




Event 2: f 

= 540 kHz 



146.9 

246.1 

4.4 

yes 

3.10 

37.8 

80.3 

16.89 

150.7 

249.5 

5.6= 

* yes 

3.20 

46.5 

77.1 

17.34 

146.9 

237.1 

5.2 

no 

3.18 

43.6 

83.5 

12.39 

150.7 

239.5 

6.4 1 

“ no 

3.26 

50.0 

79.3 

12.53 


ported a high degree of right-hand circular polarization for the 
Ulysses observations of HOM (compare their Figure 2), i.e., 
they found no evidence for significant elliptical polarization. 
Wave polarization can be easily obtained in the cold plasma 
theory of wave propagation. Stix [1962, p. 42] defines the elec- 
tromagnetic polarization as 

P = iEJE, = ( n 2 - S)/D 

where E x and E y are the wave transverse electric field compo- 
nents, n is the index of refraction, S = ( R + L)l 2 and D — 
(R — L)l 2. R and L are defined in terms of the wave real 
frequency, plasma frequency, and gyrofrequency [Stix, 1962, 
equations 19-21, p. 7]. We have calculated the polarization of 
the wave as a function of ray path. A typical example is pre- 
sented in Figure 6. Here we see that P « 1 (near linear 
polarization) near the source point, but increases rapidly as the 
wave propagates outward, toward the spacecraft. Thus the 
observation of large (>90% circular) polarization is expected 
in the cold plasma limit, consistent with the Ulysses HOM 
observations. This result is seen to follow naturally from a 
consideration of the wave characteristics along the ray path. 
Near the source region the wave frequency / ~ / RX and the 
calculated wave normal angle, tp ~ 90° as expected. As the ray 
propagates outward the index of refraction decreases, and \p 
also decreases to values «90°. Let us consider the two special 
cases, therefore, of propagation at wave normal angles near 
zero and near 90°. In the cold plasma limit it is well known for 
extraordinary waves with \p ~ 0, 

n 2 - R 

while for i// - 90°, 

n 2 ~ RLIS 

Thus, for tp — 0 

P~(R - S)/D ~ 1 

while for ip ~ 90°, 

P = ( RL/S - S)/(R - L)l 2 = (L - R)/(L + R) 

Near the HOM source region, / — / RX and ip ~~ 90°, so P < 1 
as we observe in the calculations, while near the spacecraft, i p 
~ 150° and P ~ 1 as observed. These results are not incon- 
sistent with past observations of elliptical polarization of deca- 
metric emission for which propagation elfects can be different. 


5. Discussion 

For event 1 the source locations determined from the ray 
tracing differed somewhat from those obtained in the study of 
Reiner et al. [1993a] (which assumed no refraction). For in- 
stance, at / = 540 kHz, the calculated latitude of the ray with 
the smallest L value and the latitude of the ray with the largest 
L value bracket the latitude of the source obtained by Reiner 
et al. For event 1, if the waves were launched from Ulysses 
within the estimated error limit (~4°) of the observed wave 
direction angle, the rays intercepted / RX , as required for RX 
mode emission. This was not the case for event 2 where we had 



-Y 



-Z 


Figure 5. (a) Top and (b) three-dimensional views of ray 

paths for waves launched at / - 540 kHz from the event 2 
satellite position and in the same direction relative to the 
magnetic field. In one case the 06 magnetic field model was 
used and in the other case the 04 magnetic field model was 
employed. This plot indicates differences in ray path due solely 
to the effects of the magnetic field model. 
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to modify the initial wave direction angles by more than 15° in 
order for this to occur, far outside the error estimate. 

For event 2 we investigated the effects of modifications in 
the plasma density (by comparing runs with and without a 
plasma torus) and modification of the magnetic field (compar- 
ing the 06 to the 04 model). However, in all cases it was 
necessary to modify the ray direction angles, Ak, by signifi- 
cantly more than the 4° estimated error limit in order for the 
ray to intercept / KX , a necessary condition for a source region. 

Our results indicate that refractive effects due to the Io 
plasma torus and the magnetic field are significant for HOM 
propagation. As seen in Figures 2 and 3, rays can be signifi- 
cantly refracted out of the meridian plane, indicating source 
locations that are at longitudes that differ by perhaps tens of 
degrees from the spacecraft meridian plane. Results for event 
1 (Table 2) suggest modest refraction out of the meridian 
plane from source positions located along magnetic field lines 
with L < 1 and require initial wave directions at the spacecraft 
that differ by ^3° from the observations. For event 2, however, 
we find that we must modify the launch directions at the space- 
craft by >15° (over 3 times the estimated uncertainty) in order 
tor the ray to intercept a possible source region. We suggest 
that differences between the Io plasma density and magnetic 
field models compared to actual conditions at the time of the 
Ulysses flyby may account for part of the magnitude of Ak 
obtained in our results for event 2, which occurred about 8 
hours after event 1. Temporal effects and azimuthal inhomo- 
geneities in the plasma density of the Io torus with magnitudes 
reaching a factor of 2 or more have been observed for a 
number of years fcf. Schneider et al , 1989; Schneider and 
Trauger, 1995; Desch et al., 1994; Kaiser et al 1996]. Such 
inhomogeneities have not been incorporated in our study 
which assumes azimuthal symmetry. As noted earlier, we have 
modified the entire Io torus density by factors ranging from 0 
to 5 with unsatisfactory results. Faraday rotation measure- 
ments [ Warwick and Bulk, 1964] and more recent polarization 
studies [Dulk et al ., 1994] have suggested that the electron 
density close to Jupiter may be quite low (< 1 cm 3 ). However, 
we have modified the inner plasmasphcric density ( r <4 R,) 
of the Divine-Garrett model also by factors ranging from 0 to 
5 (results not shown), but were unable to significantly decrease 
Ak for event 2. All of these modifications maintained azi- 
muthal symmetry, however. Large plasma clouds rotating 
within the torus, for instance, have not been modeled. Anom- 
alies in the magnetic field model might also explain our results 
for event 2. Connerney [1993] discusses in some detail the 
possibility of B field anomalies, especially at low altitudes 
where in situ measurements do not exist. 

These results illustrate the difficulty of accurately locating 
HOM radio sources in the Jovian magnetosphere. Even 
though in this analysis we have for the first time incorporated 
rather precise information on the incident direction of the 
radiation at the spacecraft, tracing these rays back to the 
source region requires an accurate, time-dependent and spa- 
tially inhomogeneous Jovian magnetospheric model, which is 
not realized at present. We therefore have to be content with 
an approximate HOM source location determination. 

It is significant that the largest value of L obtained in our 
ray-tracing analysis was L m;iX = 7, which is less than the L shell 
range of HOM source regions suggested by Ladreiter and Le- 
blanc [1990a] of 15 < L < 20 (auroral zone sources). The 
latter study was performed without knowledge of the observed 
wave directions that we have incorporated in our study. 


f = 940 kHz + = 153.3 /3 = 269.9 
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Figure 6. Polarization versus distance from the source ( Rj ) 
for a typical ray path. 


6. Conclusions 

We have performed a ray tracing analysis of Jovian HOM 
emission. Because of the unique direction-finding capability of 
the URAP instrument on board Ulysses, we have for the first 
time to our knowledge used measured wave directions at a 
spacecraft to initialize a ray tracing procedure. The calcula- 
tions were presented for two distinct spacecraft locations in the 
Jovian magnetosphere, where HOM emission was observed 
(events 1 and 2), and for several frequencies. The ray tracing 
calculations were performed using a model of the Io torus 
[Divine and Garrett , 1983], and two different Jovian magnetic 
field models (04 and 06). It was found that significant wave 
refraction occurred for the HOM emission due both to the Io 
torus, and the Jovian magnetic field model, however, calcu- 
lated source locations for the HOM all occur on Jovian field 
lines with 3 =£ L < 7 in agreement with the results of Reiner 
et al. [1993a]. The ray tracing results for event 2 strongly imply 
either a temporal or spatial asymmetry or anomaly in either the 
plasma density or magnetic field structure. 

Finally, calculated polarization of the emission is <5Cl near 
the source region, and —1 near the spacecraft as observed. 
This result is due, in part, to the w r ave normal angle depen- 
dence of the index of refraction. 
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